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As discovered in the quantum Hall effect, a very effective way for strongly-repulsive electrons 
to minimize their potential energy is to aquire non-zero relative angular momentum. We pursue 
this mechanism for interacting two-dimensional electrons in zero magnetic field, by employing a 
representation of the electrons as composite bosons interacting with a Chern-Simons gauge field. 
This enables us to construct a dual description in which the fundamental constituents are vortices 
in the auxiliary boson fields. The resulting formalism embraces a cornucopia of possible phases. 
Remarkably, superconductivity is a generic feature, while the Fermi liquid is not - prompting us 
to conjecture that such a state may not be possible when the interactions are sufficiently strong. 
Many aspects of our earlier discussions of the nodal liquid and spin-charge separation find surprising 
incarnations in this new framework. 
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I. INTRODUCTION 

Fermi liquid theory is the cornerstone of the modern 
theory of metals, as well as band theories of insulators 
and semiconductors. This theory - like most perturba- 
tive theories - is informed by the assumption that the ki- 
netic energy is the dominant scale. As a result, the theory 
is constructed in momentum space, where the kinetic en- 
ergy is diagonalized. This leads to strong kinematic con- 
straints which circumscribe corrections to the underlying 
free fermion behavior. In this paper, our point of depar- 
ture is a different extreme limit in which the interaction 
must be dealt with at the outset, and the kinematic con- 
straints might, consequently, be inoperative. Hence, we 
are forced to adopt a non-perturbative approach. As is 
often the case in non-perturbative problems - for exam- 
ple the quantum Hall effectfl and the one-dimensional 
electron gasa - it is advantageous to adopt a real space 
approach. Here, this also enables us to gain a vantage 
point from which to focus on the strong electron interac- 
tion. 

The past few decades have witnessed the discovery of 
a number of physical systems in which the interaction 
energy is comparable to or greater than the kinetic en- 
ergy. These materials exhibit strange behavior which is 
not readily captured within the conventional Fermi liquid 
framework. The high-T c cuprate superconductorsd are 
the most famous example, but them are certainly oth- 
ers, such as heavVifermion materialsu and high-mobility 
2DEGs at large r s U. Ironically, even 3 He - the birthplace 
of Fermi liquid theory - falls in this category!! The anal- 
ysis of such systems may require an approach of the type 



propounded in this paper. 

Following the above reasoning, we are led to search for 
a means of incorporating strong electron-electron repul- 
sion from the outset. As discovered by Laughlinjj the 
spatial separation due to the centrifugal barrier for non- 
zero angular momentum is a very effective way for parti- 
cles to lower their Coulomb energy. We consider strongly 
interacting electrons moving in the two-dimensional (2d) 
continuum, and assume that these strong interactions in- 
clude a hard-core which prevents the crossing of electron 
trajectories. Some of the resulting physics is reminiscent 
of the quantum Hall effect: pairs of particles tend to spin 
around one another. There are, however, some signifi- 
cant differences: time-reversal symmetry is not explicitly 
broken, and further, the kinetic energy is not quenched. 
Nevertheless, our investigations in the remainder of this 
paper and elsewherefl lead us to suspect that strongly re- 
pulsive electrons in the 2d continuum can form a p-wave 
(p x ± ip v ) superconductor! 

To develop a low-energy effective field theory, we first 
use only the non-crossing constraint on the fermion 
world-lines. This is a sufficieri£,condition to allow the use 
of statistical transmutationBEj to realize up and down- 
spin electrons as bosonic fields interacting with a Chern- 
Simons gauge field which attaches flux to spin. Without 
additional assumptions, we can then pass to a dual the- 
ory of vortices in the up- and down-spin bosonic fields. 
In this way, we argue .that many of our previous results 
on the nodal liquidE3E3 hold with a much wider range 
of validity. With this approach, we believe that we gain 
an unfettered view of the entire phase diagram of this 
infinitely strongly interacting fermionic system. 

Since our dual theory is of Ginzburg-Landau form, its 
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phase structure can be analyzed by considering the con- 
densation of various fields. If no vortex field condenses, 
the system is superconducting with a non-zero anau.- 
lar momentum pairing state with p x + ip y symmetryEj. 
Since vortex condensation typically implies charge or- 
dering, it is usually driven by a periodic potential or 
long-range Coulomb interactions. In their absence we 
thus conjecture that the generic state of the strongly- 
interacting system is superconducting. This is an aston- 
ishing conclusion, given the lack of a palpable 'pairing 
mechanism'. Evidently, strongly repulsive interactions 
essentially force electrons of opposite spin to "rotate" 
about one another and introduce strong superconduct- 
ing correlations. Other ordered phases result when the 
vortex fields condense. For example, a spin- or charge- 
density wave results if both the up- and down-spin vortex 
fields condense. 

A basicJeature of any superconductor is spin-charge 
separatioro. To access spin-charge separation in the 
Px + ipy superconductor within our dual Ginzburg- 
Landau formulation requires consideration of "paired" 
vortex composites. When these composite bosons con- 
dense, they can destroy the superconductivity - but spin- 
charge separation survives. Specifically, if a vortex in the 
down-spin boson field pairs with a vortex in the up-spin 
boson field and this pair condenses, translational symme- 
try is spontaneously broken by the formation of a crys- 
talline state of spinless charge e solitons. The spin sector 
is gapped except for chiral edge states, so that this phase 
is a T-violating nodal liquid (i.e. a chiral spin liquid). Al- 
ternatively, if a vortex in the down-spin boson field pairs 
with an anii-vortex in the up-spin boson field and this 
pair condenses, a transition occurs into a fully-gapped 
superconductor, such as a superconductor with tightly- 
bound pairs. If both types of vortex pairs condense, an 
analogous spin-liquid results. In each of these phases, 
time reversal invariance is spontaneously broken. 

With the inclusion of an ionic potential acting on the 
electrons, however, T-invariant phases are possible, and 
expected. Indeed, by allowing for terms in the dual 
Ginzburg-Landau theory which break rotational invari- 
ance, gapped modes in the spin sector can go soft at finite 
momentum. For a uniaxial potential, gapless modes nat- 
urally appear at two points in momentum space. At these 
two points the vortex-anti- vortex field is critical, and can 
be conveniently re-fermionized as two Dirac fields. These 
can be identified as the nodal quasiparticles of a p x su- 
perconducting phase. Similarly, an ionic potential with 
square symmetry leads to four low-energy points in mo- 
mentum space, and thereby a fourfold Dirac theory re- 
covering the spectrum of d x 2_ y 2 quasiparticles (see be- 
low and Sec. VI for a discussion of some subtleties of 
the d-wave case). Within our theory, a (very) strong 
local repulsion acting in concert with an ionic potential 
with square symmetry are the essential ingredients for 
high temperature <i-wave superconductivity. In the ab- 
sence of the ionic potential, strong T-violating pairing 
with p x + ip y symmetry is expected. This is the pairing 



symmetry in the A-phase of a superfluid 3-He filmti 

An appealing feature of our dual Ginzburg-Landau 
formulation is that it gives a clear meaning to spin- 
charge separation (and spin-charge confinement) in two- 
dimensional electron systems. Indeed, we identify an 
Ising-like Z2 symmetry which when unbroken leaves spin 
and charge separate. Spin-charge confinement is driven 
by an Ising ordering transition. 

Remarkably, although our theory is intimately tied to 
a real-space picture, Fermi surface physics is not lost, as 
evidenced by the nodal quasiparticles. The Fermi liquid 
phase itself is much more elusive! It occurs, if at all, 
only as a narrow region analogous to a (fermionic) "en- 
tangled flux liquid" in the Ginzburg-Landau theory, and 
is certainly only possible if the Z% symmetry is broken. 
However, it is unlikely that such a state can be obtained 
in the Ginzburg-Landau theory. It is here that we spec- 
ulate the relative angular momentum between electrons 
plays an important role. Indeed, (very singular) s-wave 
scattering states are known to dominate the physics of 
the Fermi liquid with hard-core radius a <C l/kp. We 
suspect that the inclusion of such singular wavefunctions 
invalidates the passage to a local and analytic vortex la- 
grangian. See Sec. VI for further discussion of this point. 

The phase diagram which results from this analysis 
contains a plethora of fascinating states, including su- 
perconducting states of pairing symmetry p x , p x + ip y , 
d x 2_ y 2 1 d xy , and their quantum disordered counterparts. 
These states are characterized by a separation between 
the characteristic scales of the charge and spin and even 
- when the Z2 symmetry is unbroken - true spin-charge 
separation. A particularly salient aspect of this phase 
diagram is the miniscule domain, and we speculate per- 
haps even absence, of the Fermi lipoid state. This seems 
to support the claims of AndersonEj and co-workers that 
Fermi liquid theory breaks down in strongly-interacting 
two-dimensional systems (though it does not support 
their bold claims that Fermi liquid theory breaks down 
even at weak-coupling) . The proofs p£ the stability of 
Fermi liquid theory at weak-couplingtll presumably do 
not apply to the strongly-interacting limit which we have 
considered. 

In section II, we first discuss the statistical transmuta- 
tion which obviates the need for a local pair field. Con- 
structing a dual theory, we describe the standard antifer- 
romagnetic and charge-density-wave states which result 
from the condensation of single vortices. In section III, 
we discuss paired vortices and the Z2 symmetry which 
distinguishes their condensation from that of individual 
vortices. When this symmetry is unbroken, spin and 
charge separate. In section IV, we discuss the phase 
diagram which results from the condensation of paired 
vortices. This phase diagram revolves about a p x ± ip y 
superconducting state. In section V, we show how time- 
reversal (T) invariant superconducting states such as p x 
or d xy can arise in this model. The momentum-space 
structure and concomitant phenomenology of BCS-like 
d xy or d x 2_ y 2 superconductors is recovered. We find 
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that the physics of the nodal liquid reappears in a new 
guise: the nodons are vortices in a vortex field whose 
fermionic statistics result from their interaction with a 
Chern-Simons gauge field. In chapter VI, we arrive, ulti- 
mately, at a phase diagram which is the synthesis of ideas 
of duality and vortex condensation common to field theo- 
ries of the quantum Hall effect as well as our earlier work 
on nodal liquids, but is almost entirely orthogonal to the 
underlying conceit of Fermi liquid theory. 

II. FERMIONS, FLUX ATTACHMENT, AND 
DUALITY 

We focus throughout on spinful electrons moving in 
the two-dimensional continuum, interacting via a spin- 
independnet interaction. We assume that the electron- 
electron repulsion is strong enough that no two electrons 
can ever be coincident. Precisely this "hard-core" con- 
straint makes it possible to transform the interacting two- 
dimensional electron gas into a mathematically equiv- 
alent system of interacting spinful bosons, by attach- 
ing "statistical" flux with an appropriate Chern-Simons 
gauge field. Such a "bosonization" scheme for 2d spinless 
electrons has been particularly illuminating, in the con- 
text of the fractional quantum Hall effectEJiMi. With spin 
there is considerable freedom in how one attaches the 
flux tubes to convert fermions into bosons. We adopt a 
scheme in which flux is attached to the spin of the elec- 
trons, and define 

c a (r) = b a (r) exp[ze Q f 9(r - r')2S z (r% (1) 

Jr' 

with a "charge" e-f = 1 and ej_ = — 1. Here 

S*(r) = [nr(r)-n;(r)]/2 (2) 

is the z-component of the spin density operator with 
n a = cl t c a = b\ t b a (no sum on a), and 0(r) denotes 
the angle that r makes with the rr-axis. The boson oper- 
ators satisfy canonical commutators, [b a (r), b^(r')] = 
for r/r'. Due to the non-crossing constraint, the "on- 
site" commutators need not be specified. 

An advantage of the above scheme for flux attachment 
is that with zero total spin, Sf ot = (as assumed here- 
after), the statistical flux "seen" by the Chern-Simons 
bosons vanishes on average. For spinless electrons this 
happy situation requires the presence of a strong exter- 
nal magnetic field (as in the FQHE, c.f. Ref. As 
for Abelian bosonization in one spatial dimension^, the 
above choice of a spin quantization axis masks the un- 
derlying spin-rotational invariance. But as we shall see, 
it is possible to restore explicit SU(2) symmetry by a 
subsequent "refermionization" . 

After transforming to boson operators, the partition 
function for the Hamiltonian of 2d interacting electrons 
can be expressed as a functional integral over bosonic 



fields and a statistical gauge field, a M , with associated 
(Euclidean) Lagrangian density: 

£ = K(d T - ie a a°)b a - ^—b* a (y - ie a a) 2 b a - C cs . 

(3) 

In the "Coulomb gauge" diet 1 = the Chern-Simons term 
is C cs — (i/2ir)a (eijdtaj), but can be cast into a more 
convenient gauge invariant form: 

£cs{au) = i-^£^u\a^d v a\. (4) 

The form of the electron interaction term (not shown) is 
unchanged under "bosonization" due to the equivalence 
of the fermion and boson densities: c^c = wb. 

We now implement, the standard 2+1-dimensional du- 
ality transformationlMj which exchanges bosons (the 6 Q 's) 
for vortices in the bosonic fields - arriving at a description 
in terms of vortex field operators, denoted $ Q . To illumi- 
nate this, it is instructive to briefly consider an alternate 
representation in terms of boson world lines: 

£*>i = \{J C Z) 2 + J>» - £c S K). (5) 

Here J" denotes a bosonic three-current (with /i running 
over 2+1 space-time coordinates) for spin component a. 
The first term measures the length of the space-time 
world-lines and represents the kinetic energy. The Chern- 
Simons coupling generates a sign change when two world 
lines exchange, transforming to fermions. To implement 
duality, these three-currents are expressed in terms of two 
gauge fields, <z^, one for each spin component: 

- ^ vX d v a a x . (6) 

Z7T 

In this way, charge conservation (d^J" = 0) is automat- 
ically satisfied. The dual Lagrangian can be obtained 
by inserting this expression into C w h an d integrating out 
the Chern-Simons field, a M . Electron charge quantization 
is implemented by the vortex operators, which are 
minimally coupled to a™. 

The final dual form consists of two Ginzburg-Landau 
theories, coupled via a Chern-Simons term: 

Cdual = C GL(<$> a , <) + Ccs(al - oji), (7) 



CGrt*,aJ = - m M )$| 2 + V(*) + \{U„) 2 . (8) 

with a "potential" that can be expanded as V($) = 
r|$| 2 +w|$| 4 + .... 

Since the Chern-Simons term only involves spin cur- 
rents, it is extremely convenient to introduce charge and 
spin gauge fields: 
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< = < + «i t ; = (9) 

(and corresponding field strengths, As with 

Abelian bosonization in one dimension, charge and spin 
currents defined by J? = + and Jg = J.J — J^, 
are given by derivatives of the charge and spin fields, 
respectively: 

J £ = ^-^uxd u a p x ; J° = ^-e^\d u a a x . (10) 

Z7T Z7T 

Longer range Coulomb interactions can be readily incor- 
porated by adding a term bilinear in the electron charge 
density: eijdiCLj. In this dual representation, $t creates 
a vortex in the electron wavefunction - effectively increas- 
ing the angular momentum of all spin a electrons by one 
unit. An important feature of Cduai is that there are pre- 
cisely as many positive as negative circulation (electron) 
vortices (for each spin species), implying a relativistic 
form for $ Q . 

It is instructive to briefly mention how these dual fields 
couple to an external electromagnetic field, A^. As usual 
couples directly to the total electrical current, so 
that from Eq. [h] one has: 

C A = — A^AC^a?. (11) 

It is also convenient to introduce an external "spin" gauge 
field, A" , which couples to the total (z-component) spin 
current, J£: 

£a = (12) 

The dual Ginzburg-Landau representation can be 
fruitfully employed to describe various possible phases 
of spinful 2d electrons satisfying the 'non-crossing' con- 
straint. For instance, imagine a phase in which the (elec- 
tron) vortices are absent in the ground state (except as 
virtual fluctuations), which corresponds to taking r large 
and positive in the above Ginzburg-Landau description. 
Being massive, $ Q can be safely integrated out, leaving 
an effective theory: 

£eff = l(fp 2 +£ cs (a;), (13) 

which describes massless charge fluctuation, and spin 
fluctuations gapped (in the bulk) by the Chern-Simons 
term. This is a superconducting phase, which can be ver- 
ified by noting that the pair field operator (efCjJ creates 
a (2+1 space-time) monopole of strength two in the field 
= e iiV \cl v a p x . Since the gauge field a£ is massless, the 
energy cost to make a free monopole is finite (monopoles 
interact via 2+1 Coulomb force), so that the pair-field 
exhibits true ODLRO. Due to the Chern-Simons term, 
this superconducting phase exhibits a quantized Hall 
"spin-conductance", a s xy = 1, a signature of a sponta- 
neous breakdown of time reversal invariance. This fol- 
lows by noting that the Lagrangian with spin "gauge" 



field, C e ff + £% depends quadratically on both a£ and 
a^, so that they can be integrated out to give, 

C eff +C A = -ia s xy -^e^ x A°d v AZ, (14) 

with a xy = 1. Following the analysis in Ref. |l^, one 
can readily verify that a 2d BCS superconductor with 
Px + ipy pairing symmetry has precisely such a value for 
the quantized spin conductance (also see below in Sec- 
tion IV) . The spin state of the pair is then presumably a 
triplet with s z = 0. This is the phase of a 2d superfluid 
3-He A film. 

Before discussing spin-charge separation, which is a 
generic property of a 2d superconductor, it is instructive 
to consider phases described in the dual theory when vor- 
tices created by $ Q proliferate, and condense: ($ Q ) 0- 
As we shall see, in contrast to the p x + ip y superconduct- 
ing phase, these phases typically exhibit crystalline or- 
der, spontaneously breaking translational symmetry. To 
see this, note that upon vortex condensation, the dual 
"flux" Cijdia" is quantized in units of 2ir, which corre- 
sponds to quantization of charge in units of the electron 
charge e. By analogy with the Abrikosov flux lattice 
phase of a Type II superconductor, one expects a break- 
down of translational symmetry with spin up (and down) 
electrons forming an ordered lattice. Depending on the 
relative phase between the density wave of spin up and 
down electrons, this will be either a charge density wave 
(CDW) state or an antiferromagnet (AF). In the pres- 
ence of a commensurate background periodic potential 
from the ions in the solid, one expects these density wave 
states to lock, resulting in insulating behavior. 

These "crystalline" phases can presumably be energet- 
ically stabilized by a longer-range repulsive interaction 
between the electrons, in addition to the 'non-crossing' 
constraint (which is required to make our 2d bosonization 
scheme legitimate). In their absence, our dual Ginzburg- 
Landau representation, (fortified by subsequent analysis 
below and physical reasoning) strongly suggests that the 
predominant ground state is the p x + ip y (or p x — ip y ) 
superconductor. Given that spin-up and spin-down elec- 
trons prefer a state of non- vanishing relative angular mo- 
mentum to minimize Coulomb repulsion, an £ — ±1 or- 
bital angular momentum state is clearly favored by the 
kinetic energy. As such, it seems that incorporating local 
Coulomb repulsion by forcing electrons of opposite spin 
into a relative angular momentum state is a very effective 
electronic mechanism for high temperature superconduc- 
tivity! 

III. SPIN-CHARGE SEPARATION 

The phases described above are not the only possible 
phases for 2d electrons satisfying the 'non-crossing' con- 
straint. Composite order parameters can also condense, 
thereby leading to charge- and/or spin- insulators. We 
focus on the combinations 
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$„ = * T $ i; 4> ff = $ T ${, (15) 

which, as we shall see, are exceedingly interesting from 
a phenomcnological standpoint. These order parameters 
can condense without breaking the Zi symmetry 

$ u _>_$ Uj (16) 

so we can have (& P ,a) ^ while <&|,j = 0. 

In the following we presume that $ p and <& a describe 
the soft modes at low energies, and that $ Q remains mas- 
sive. As we shall see, this leads naturally to a separation 
of low energy spin and charge degrees of freedom. Our 
motivation for this is two-fold. Firstly, spin-charge sepa- 
ration is a generic property of a superconductor such as 
the p x + ip y state discussed above, and it is instructive 
to exhibit this separation within the present Ginzburg- 
Landau framework. But secondly, in many Mott insula- 
tors of interest the charge degrees of freedom freeze out at 
much higher energy scales than the energies on which lo- 
cal moments and spin order develops. This is typified by 
the undoped cuprate materials, with insulating behavior 
setting in on the scale of electron volts (the "Hubbard" 
U) much higher than the antiferromagnetic ordering tem- 
perature. In order to capture these two very different en- 
ergy scales within the present framework, it is essential 
to transform to the charge and spin vortex fields, <& p and 
<f> CT . Indeed, in the description of the antiferromagnetic 
insulator discussed above driven by condensation of <i> Q , 
charge ordering and local moment formation necessarily 
take place on the same energy scale, since the dual flux 
tubes in these vortex fields are electrons carrying both 
charge and spin. 

Under the assumption that both fields $ Q remain mas- 
sive, one can write down an effective theory for the soft 
modes $ p<a by integrating out ^t.j.- Below we illustrate 
how this can be done, by regularizing the theory on a 
lattice. But more generally, the form of the effective the- 
ory is essentially dictated by symmetries, involving three 
contributions: 

£eff=£p + £a+£int, (17) 

with a charge sector, 

Cp = \\(P»- 0<!> P \ 2 + r P \^ P \ 2 + u p \$ P \ 4 

+ \{fP 2 + ^A M d v a% (18) 

a spin sector, 

Ha = = \\{d P - ia;)$>a\ 2 + r CT |$ CT | 2 + u a \<5> a \ A 

+\{f%) 2 +i^uxald v al (19) 

and sub-dominant interaction terms involving many 
derivatives (see eg. below). The charge sector has the 
Ginzburg-Landau form, with minimal coupling to the 



charge gauge field, a£, and A p is the physical electro- 
magnetic potential. The Chern-Simons term lives solely 
in the spin sector. 

Some insight into the genesis of such a Lagrangian may 
be obtained by considering a lattice version of (ph and 
dropping the for simplicity. Writing <I> Q = e a , we 
have 

S=5>os(0?-0?), (20) 

(i,3) 

where i and j denote sites of a 2d (say square) lattice and 
a sum over a is understood. We now introduce charge 
and spin fields p ' a : 

u = 5(0'±<n + |» (2i) 

where s = ±1 is an Ising "spin" variable. By introducing 
s, we can treat 8 p,cr as angular variables since the action 
is invariant under d p,rT — > 9 p - a + 2tt, s — > — s. The action 
can then be rewritten as: 

S = 5> Sj cos 1(0? - %) cos l(0f - 0f), (22) 

since sin^Sj — Sj) = and cos^(si — Sj) — SiSj. Let 
us now consider the effect of integrating out the Sj's. If 
we are in the symmetric phase in which the Zi is un- 
broken, this can be done perturbatively, as in the high- 
temperature expansion for the Ising model. 

Z = J2 e^ fiJsiS} 

s,=±l 

Si=±l n ■ \ ij J 

(23) 

To leading order in /?, which corresponds to decoupled 
free spins, one has (siSj) — Sij, which implies an effective 
action of the form: 

Seff = £ (1 + COs(0f - %)) (1 + COB(0? - Of)) 

(i,3) 

= £cos(0?-0?) + cos(0f-0J) 

(id) 

+ cos(0£-^)cos(0f-0J). (24) 

Upon making the identifications $ p = e l6lp , & a = e l6a 
and restoring the gauge fields (minimally coupled) the 
first two terms are seen to be lattice versions of the con- 
tinuum Ginzburg-Landau theories in C p and respec- 
tively. The last term generates a gradient interaction 
term between the charge and spin sectors. 

When the Zi symmetry is unbroken, as it is by as- 
sumption in (|l7|), spin and charge separate, as we now 
argue. In the p x + ip y superconducting state, which can 
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be described by either (||) with r-f^ > or fll7| ) with 
r p ,a > 0, the low-energy excitations are the gapless su- 
perfluid mode, a?, which carries charge but no spin. At 
finite energy, there are also the quanta of <$> a , which are 
fcrmionic and carry spin-1/2 by virtue of their coupling 
to the Chern-Simons gauge field, a a . They do not couple 
directly to the electromagnetic field, so we assign them 
quantum numbers q — 0, s = 1/2. As we shall see in 
the next section, these neutral fermionic spin 1/2 excita- 
tions are the p-wave analog of "nodons", introduced in 
reference |ll| for a d-wave superconductor. 

When r p < and $ p condenses, the dual flux e^did? 
becomes quantized into "flux-tubes", by direct analogy 
with a Type II Ginzburg-Landau superconductor. Each 
one of these dual 'flux tubes" carries one unit of electric 
charge, but no spin. _We refer to these spinless charge 
e solitons as "holons"lEl. We thus see that provided the 
Z2 symmetry in Eq. |l6] is unbroken, spin and charge are 
separated. On the other hand, when the Z 2 symmetry is 
broken by the condensation of $f,j., the spin and charge 
are confined. Since (fc^ couples to a£ ± this conden- 
sation locks the spin and charge together, leaving only 
the electron in the spectrum. 

To summarize, states of higher symmetry have less re- 
stricted spectra. The original dual representation in Eq. 
|§| has a U(l) x U(l) gauge symmetry, corresponding to 
independent rotations of <J>j and 

$ Q -» $ a e <A °; 0° -» al + dfj,A a , (25) 

with two arbitrary functions A a . This gauge symmetry 
emerges when the conserved electron three-currents are 
expressed as a curl of the gauge fields, a™. Breaking 
down this large symmetry corresponds to "localization" 
or "quantization" of charge and/or spin. When the full 
symmetry is completely broken, both charge and spin 
become quantized together, and all the excitations have 
quantum numbers of the electron, with | + s an inte- 
ger (as in the antiferromagnetic insulator mentioned in 
Section II) . But if this symmetry is only partially broken 
by condensation of <E> P and leaving an unbroken Z2, 
both charge (e) and spin (1/2) become quantized, but 
excitations exist with any combination of these quantum 
numbers. This will be nicely illustrated in the next Sec- 
tion where we employ the dual description, Eq. to 
examine the properties of some T-breaking spin-charge 
separated phases. 

IV. PHASES WITH BROKEN T 

Having established the form of Eq. [ll] under the as- 
sumption of an unbroken Z2 symmetry, and the concomi- 
tant spin-charge separation, we explore possible phases 
which emerge from this effective theory. We first focus 
on phases which (spontaneously) break time reversal in- 
variance. As we shall see, these emerge naturally for 



electrons in jellium, moving in the absence of ionic po- 
tentials. In the following sections we consider the effects 
of ionic potentials which break rotational invariance and 
naturally drive transitions into time reversal invariant 
phases. 

As already discussed, the phase in the absence of ei- 
ther spin or charge vortices is a p x + ip y superconductor. 
Since the charge and spin sectors have effectively decou- 
pled under the assumption of the unbroken Zi symmetry, 
it is possible to consider them separately. If the vortices 
in the charge sector proliferate and condense, (<3> p ) 7^ 0, 
the field e^diUj becomes quantized in "dual" flux tubes, 
as discussed above. Each of these "flux" tubes carries 
charge e, but no spin. These charge e spinless "holons" 
are expected to crystallize, by direct analogy with the 
Abrikosov flux-lattice. This charge e crystal will pre- 
sumably lock to any underlying ionic potential. With 
gapless spin-carrying edge states still present, this electri- 
cally insulating phase is the p x + ip y analog of the nodal 
liquid. Once again, energetic stabilization of this crys- 
talline phase presumably requires the presence of appre- 
ciable Coulomb repulsion between electrons on the scale 
of the mean electron spacing (in addition to the 'non- 
crossing' constraint). 

But suppose the spin vortices condense, in the absence 
of charge vortices? Due to the dual Anderson-Higgs 
mechanism, becomes massive rendering the Chern- 
Simons term ineffective, and leading to a spin-gap both in 
the bulk and at the edge. This implies "spin-insulating" 
behavior with <j s xy = 0. How can we understand this 
fully spin-gapped superconducting phase? To this end it 
is convenient to briefly consider a BCS description of the 
quasiparticles in a p x + ip y superconductor: 

Hbcs = e kc{ a c ka + A fc cj cT c t _ fci + h.c, (26) 

k 

with dispersion e& = (k 2 /2m e ) — [i and gap function 
Afc = V/\(k x + iky). In terms of a two-component spinor, 
-01 (r) = C|(r) and ^(f) = c[( r ), this can be rewritten 
in the form of a Dirac equation with (Euclidean) La- 
grangian: 

Cbcs = ^[d T + r z ((-d 2 /2m e ) — fi) + iv A T j dj]iP, (27) 

with j = x,y. This gives the usual BCS quasiparticle 
dispersion: = ±y 'e k + {v&k) 2 . As in Ref. [ll| one 
can define a gauge invariant charge neutral quasiparticle 
(a "nodon"), by transforming ip — > exp(ir z ip/2)ip, with ip 
the phase of the complex pair-field. Spin and charge are 
thereby separated, with the z-component of spin being 
the conserved U(l) "charge" in the Dirac theory: 2S Z = 
ip^tp. Since the source field A° couples to the conserved 
spin current, it can be readily incorporated into the above 
Dirac equation via a "minimal coupling" prescription: 

d, - - iA°. 

In the presence of a boundary, say at y — with 
boundary conditions c a (x,y = 0) = 0, one can readily 



G 



show from the above Dirac theory that a chiral fermion 
edge mode exists only for positive chemical potential, 
fi > 0. In this BCS limit one clearly has a xy = 1. But 
at very strong coupling when /i changes-sign, the ground 
state changes to a paired "molecular" E3 limit with zero 
a xy = 0. Right at the transition, there are gapless bulk 
quasiparticle excitations described by a massless Dirac 
theory (at /i = in Eq. ^7|) with a "node" at zero mo- 
mentum. To access the molecular limit of the p x + ip y 
superconductor presumably requires a very strong (and 
unphysical) attractive interaction between electrons, en- 
abling up and down spin electrons to form a finite angu- 
lar momentum bound state (with I = 1). The attractive 
interaction must overcome the centrifugal repulsion be- 
tween the two electrons (present due to the 'non-crossing' 
constraint). 

A direct connection between the molecular p x + ip y su- 
perconductor and the phase described by the dual theory 
when the spin vortex condenses, (3><j) 7^ 0, can be estab- 
lished by re-fermionizing the spin sector of the Ginzburg- 
Landau theory and showing its equivalence to the Dirac 
theory Eq. To illustrate this we instead bosonize the 
Dirac theory. In the BCS limit of the p x + ip y super- 
conductor with fi positive, the massive relativistic Dirac 
fermion, ip, can be converted to a relativistic boson, $, 
via a Chern-Simons transformation: 



C 



BCS 



+£ 



ia u 



1A 



{AD- 



+ A/ 2 |$| 



urn 



(28) 



Here M > can be equated with the Dirac mass /i. 
Indeed, the spectrum of this massive boson field, <&, is 
ujf, = ±v M 2 + k 2 - the same form as the BCS quasipar- 
ticle dispersion, E^. Since $ is massive it can be safely 
integrated out. The only remaining dependence on the 
source field A° is through the last Chern-Simons term, 
which has been included to give the correct result for the 
spin Hall conductivity: a xy = 1. 

Remarkably, Eq. ^8] precisely coincides with the spin 
sector of the dual Ginzburg-Landau theory, C a in Eq. IE . 
Indeed, with inclusion of the source term, C a A in Eq. 12, 
the full Lagrangian in the spin sector, £(^4„) = C a + C 
can be conveniently rewritten by shifting — > a° L 
as, 



A 



iA^ a \ 2 +r a \^ a 



with r a negative. This corresponds to taking /1 negative 
and entering the molecular limits of the p x + ip y super- 
conductor. The critical point at r a = 0, with massless 
but uncondensed is equivalent to the single massless 
Dirac field (with \i = 0) centered at zero momentum. 
Without recourse to re-fermionization, vortex condensa- 
tion =/= directly implies a mass for and a van- 
ishing spin Hall conductivity, a xy = - the correct value 
for the p x + ip y molecular superconductor. This internal 
consistency gives us some confidence in the more general 
validity of the dual Ginzburg-Landau formulation. 

An alternate route from the BCS to molecular limit is 
possible by implementing a duality transformation on the 
bosonic theory C ai which interchanges the two phases. 
This can be achieved by expressing the bosonic three- 
current for the conserved spin ($t $ CT = ip^-ip) as the curl 
of a gauge field, a M , and integrating out a^. After shifting 
a M — > ctfj + A"^ one thereby obtains, 

C Dua i = ^\{dn- ia„ - iA^)4>\ 2 + r^ 2 + u|0| 4 + £ cs (a M ). 

(30) 

Here cj> creates a vortex in the field Notice that the 
dual theory has the same form as C a in Eq. [2^, except 
for the absence of the Chern-Simons term in A a . Under 
duality, the "ordered" phase with r CT < (and 7^ 0) 
maps into the "disordered" phase for <j) with r^ > 0. In 
this phase (the molecular limit) the dual theory correctly 
predicts a xy = 0, due to the absence of the A° Chern- 
Simons term. In terms of the original Dirac field this 
duality is a particle/hole transformation, tjj — > ijv , which 
changes the sign of the Dirac mass, /1 — * — fi. The self- 
dual point where both <f> and $ CT are critical, corresponds 
to the massless Dirac theory. 

From the p x + ip y molecular superconducting phase 
with o~ xy = 0, it is possible to also proliferate and con- 
dense the charge vortex: ($ p ) 7^ 0, which describes a 
fully spin-gapped crystalline phase of spinless charge e 



(29) 



This is identical to Cbcs under the identification: <I> — 
^cr, a M = a° and M 2 = [i 2 = r a . The upshot is that a 
simple re-fermionization of the spin sector, C a , gives di- 
rectly the BCS quasiparticle Lagrangian Cbcs m Eq. |2(]. 
The spin carrying but charge neutral vortex field, is 
thus seen to be equivalent to a "nodon" destruction op- 
erator. 

By such a re-fermionization procedure, we can infer 
the properties of the vortex condensed phase, ($<j) 7^ 0, 
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Fig. 1: The phase diagram in the r p — r a plane. Here, 
nodal liquid is denoted as NL, and a xy is the spin Hall 
conductivity. 
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In all four of the above phases (depicted schematically 
in the r p -r a plane in Fig. 1) time reversal invariance is 
spontaneously broken, reflecting the underlying p x + ip y 
"pairing" symmetry. In each case, this symmetry break- 
ing is taking place in the spin sector of the theory. The 
charge sector described by the simple Ginzburg-Landau 
theory (with no Chern-Simons term) is manifestly time 
reversal invariant. A natural question that arises is 
whether the spin-charge separated dual vortex theory in 
Eq. [l?] can describe time reversal invariant phases, such 
as a d x 2_ y 2 superconductor or nodal liquid. Clearly a 
d x 2_ y 2 phase requires the breakdown of rotational invari- 
ance, either spontaneously or by the presence of an ionic 
potential. Moreover, since the d x 2 _ y 2 superconductor ex- 
hibits gapless nodal excitations with finite momentum 
(on the Fermi surface in weak coupling BCS) , it is neces- 
sary to access Fermi surface physics at non-zero momenta 
in the dual formulation. This can be readily achieved as 
we now discuss. 



V. "FERMI SURFACE", P x , AND D X 2_ Y 2 
PHASES 

A. p x state 

In this section we imagine introducing an ionic po- 
tential with the lattice symmetry. The simplest case is 
a uniaxial perturbation, as might be appropriate, e.g. 
in a quasi-one-dimensional superconductor such as the 
cuprate ladder materials (the more interesting general- 
ization to square symmetry will be returned to later). 
For simplicity, we take the symmetry axes along x and 
y. Physically, it is clear that such a potential favors the 
formation of a real (non T-breaking) paired state such 
as p x or p y . This can be seen by considering the two- 
body problem deep in the molecular limit. In the pres- 
ence of the lattice potential, angular momentum is no 
longer a good quantum number (only discrete ±7r rota- 
tions and reflections are symmetry operations) and the 
p x ±ip y states will generally be mixed. Coupling the two 
via a small "tunneling" perturbation, as appropriate for 
a weak ionic potential, splits the two initially degener- 
ate levels into non-degenerate p x or p y eigenstates. The 
system then condenses into the lower of the two states. 
In the BCS limit the effects are more subtle, as we now 
illustrate. 

As alluded to in the previous section, a distinguish- 
ing characteristic of the BCS theory of unconventional 
superconductors is the presence of gapless quasiparticle 
excitations at the intersections of the nodal lines of the 
pair wavefunction with the Fermi surface. For the p x 
case, this occurs at two antipodal points in momentum 
space k — ±Ky, and the resulting gapless quasiparticles 
can be cast into the form of two Dirac species (see be- 
low) . The multiplicity of Dirac fermions signals the emer- 
gence of two conserved U(l) charges in the low-energy 



limit. Essentially, because of phase space restrictions, 
the spin at each nodal point is separately conserved - 
equivalently, one may view the two conserved charges as 
spin and (quasi-)momentum, the latter obtaining due to 
the absence of momcntum-nonconserving umklapp pro- 
cesses for generically situated nodal wavevectors. 

In the dual theory of Eq. 17, total spin conservation 



is manifest owing to the transverse nature of the spin 
^uyxdvoR. Momentum conservation is 



current J° 

H 1 



also clear, but apparently on a very different footing - 
it follows from the space-time Galilean invariance of the 
model. Apparently, to describe the BCS limit, a con- 
nection must be made between the internal (U(l) or 
SU(2)) symmetry of the model and the external space- 
time translational symmetry. It is quite remarkable that 
such a connection can indeed be clarified, as we now 
show. 

To see how the finite momentum physics can emerge 
in the dual theory, first recall the quasiparticle structure 
of the p x + iepy superconductor in the BCS limit. This 
is described by Eq. with /i > and Ak = v\k x + iek y 
(only the symmetry, and not the particular form of this 
gap function is important in what follows). The low en- 
ergy quasiparticles occur near k x = 0, k y = ±K, where 
K = y/'2m e fJb. To focus on these two regions in mo- 
mentum space, next define "slowly- varying" quasiparticle 
fields f sa via 



c t( x ) ~ ~^ (-^T - if *l) 



JsKy 



s=± 



ft j e isK v, 



One thereby obtains the Lagrange density 



C 



p+iep — / J 
s=± 



fl [d T + ivAT x d x + isv F T v d y 



(31) 
(32) 

(33) 

where vf = K/m e , m — eK and, following Eq. we 
have introduced Pauli matrices f which act in the a =| , J, 
space. Eq. |3^ has the form of two decoupled Dirac equa- 
tions, and thereby displays two manifest U(l) symme- 
tries. Note from Eq. |3^ that the f [ fields are defined as 
hermitian conjugates relative to the /j fields in Eq. pH) . 
This implies that the overall U(l) rotation, f s — ■> e lx f s , 
corresponds to spin (S z ) conservation. The relative U(l), 
f s — > e lsx f s , embodies instead translational symmetry, or 
momentum conservation. 

The conserved densities are determined by Noethcrs 
theorem: 



S- — f-f- — ^-KX C -K\ 



(34) 
(35) 



These are chiral spin densities, closely analogous to the 
chiral densities encountered in one-dimensional Fermi 
systems. The total spin density S z = + S^. 
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We are now in a position to bosonize the p x + iep y 
BCS model. The two flavors of Dirac particles f± may 
be traded for two complex boson fields <I>± by attaching 
flux in a variety of ways. The most natural, however, 
is to introduce a single gauge field, attaching one flux 
quantum to each overall (S z ) U(l) charge. This choice 
ensures both that all the resulting bosonic fields commute 
at different space-time points and that the gauge field has 
the same physical meaning as the defined previously. 
In particular, 



S z = T^V X (f. 

2vr 

With this choice, the bosonized Lagrangian is 



(36) 



C p+lep = |A)$ S | 2 + «| l^l^l 2 + V% \D 2 <S> S \ 2 

S 

+W(<S>±)+£ cs {a°)-£ cs (A' T ), (37) 

where = — i(a° + A^). Here we have included an 
external spin gauge field A° for bookkeeping purposes. 
The potential W is dictated by symmetry to take the 
form 



W(#±) =m 2 [|$+| 2 + |$_| 2 ] 



[|$ + | 2 + |<&-| 2 ] 2 + «|<MM 



(38) 



The final term in Eq. |37] is dictated by requiring a s xy = 1 
in the p x + iep y phase where m 2 > 0. 

Remarkably, the selection of two such non-zero "nodal" 
points occurs fairly naturally in our dual Ginzburg- 
Landau theory once uniaxial (rectangular) anisotropy is 
included. Consider the modified version of Eq. 

£ a = |A>^| 2 + |£>i$<x| 2 + c\D 2 $ a \ 2 
d 



(39) 



For simplicity, we have included only a single symmetry- 
breaking term, the coefficient c < I, which favors fluctu- 
ations along the y axis over the x axis. The coefficient 
d > is included for stability purposes. 

As c is decreased from one (zero anisotropy), the en- 
ergy cost for fluctuations of with spatial variations 
along y becomes more and more reduced. When c 
changes sign and becomes negative, the lowest energy 
fluctuations bifurcate away from the origin in momen- 
tum space and move to two points k = ±Ky, with 
K = y/\c\/d. From this point on, it is appropriate to 
focus on the low-energy field configurations, viz 



$_e 



-iKy 



(40) 



The physical meaning of Eq. ^ is clear from the above 
"reverse engineering" of the field content of the p x + iep y 
superconductor - compare with Eqs. |3l]-[32| Naively in- 
serting Eq. 40 in Eq. p9 and neglecting rapidly oscillating 



terms in the usual way gives an effective Lagrangian for 
the $> cr ± fields. This has precisely the form of Eq. pTL 



with v\ = 1 + |c|, Vp = 2|c|, m 2 = r a — c 2 /2d, u = u a , 
and v — 2u a . These values (particularly u and v) should, 
however, not be taken too seriously, as they certainly de- 
pend upon the simplistic treatment of fluctuations and 
higher-order terms. 

With this identification in hand, we conclude that 
Eq. provides a unified description of the p x ±ip y and p x 
states in an intrinsically anisotropic system. For c < 0, 
the equivalence to Eq. |37] allows a refermionization to the 
form in Eq. |33|. It is natural to associate the critical point 
of these equations (m = 0) with the spin structure of the 
p x superconducting state and its nodal- liquid/holon lat- 
tice counterpart. The refcrmionized double Dirac form in 
Eq. 33 is our primary result for the uniaxially anisotropic 



model. 

Issues of time-reversal symmetry merit some discus- 
sion. Ideally, a general formulation of the problem should 
contain T-non-invariant terms only through spontaneous 
symmetry breaking. However, generically both Eq. |37| 
and Eq. |3^ break T. In the fermionic formulation, Eq. |33j , 
fortunately, the explicit symmetry breaking can be eas- 
ily restored by requiring m 2 — 0. Furthermore, all T- 
preserving perturbations of this form can be showja to be 
irrelevant, so that the p x state is locally stable.! From 
the point of view of Eq. [?7|, this is remarkable. Indeed, 
a direct mean-field analysis would suggest that a gapless 
state occurs only along the critical line m 2 — 0, requir- 
ing tuning of a parameter. These conclusions can be 
reconciled by noting that constraints upon Eq. [37] in a 
T-invariant system are not at all obvious. Only because 
of the ability to refermionize are we able to identify the 
critical line m 2 = as containing a T-invariant manifold. 
If time-reversal is explicitly or spontaneously broken, a 
critical state is indeed non-generic, and the spectrum of 
Eq. |37| correctly reproduces that of the Dirac theory for 
small non-zero m 2 > (i.e. p x +iep y with e <C 1). Similar 
subtleties render the analysis of Eq. [37] problematic for 
m 2 < 0. We suspect that this region represents rather 
more exotic T- violating states, and do not consider it 
further. 

Another physical route away from the BCS-p^ phase 
is via a molecular p x state. Indeed, as we have argued 
above, deep into the molecular limit, uniaxial anisotropy 
guarantees a p x state. To tune through such a transi- 
tion, we may imagine introducing a finite-range attrac- 
tion into the 'non-crossing' model which favors tighter 
pair binding. In such a model the p^-BCS state natu- 
rally undergoes a transition into a molecular p x state as 
the attraction is increased. From the conventional BCS 
point of view, we would expect such a transition to be de- 
scribed by taking the chemical potential through zero in a 
p x quasiparticle Hamiltonian, i.e. Eq. with = vk x . 
For this model, as fj, passes from positive to negative, 
the quasiparticle nodes converge and coalesce at the ori- 
gin, becoming massive for // < 0. Precisely at the critical 



point, one expects a spectrum lu 



X/ Ami 



This 
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unconventional non-Lorentz invariant form is precisely 
what is obtained at small wavevectors from Eq. [39| at the 
point where c = r a = 0. Just as the gapless p x state oc- 
curs only along a line in mean-field theory but comprises 
a phase in the gauge model, it appears that the critical 
state at c = r a = (see Fig. 2) in reality forms a phase 
boundary despite appearing as a multicritical point in 
the mean- field treatment. 

Given the existence of the p x molecular state, there 
must be another critical line separating this from the 
p x ± ip y molecular phase. As neither state contains 
gapless spin excitations or possesses a non-zero a xy , we 
have been unable to discriminate between them within 
the GL model. Viewed from the point of view of self- 
consistent BCS theory, the transition would appear to be 
a simple first-order level crossing of two-particle bound 
states. Some preliminary modelingp suggests that a con- 
tinuous transitions in the universality class of the quan- 
tum transverse-field Ising model is also possible, the Ising 
order parameter reflecting the presence or absence of T- 
breaking. The full proposed phase-diagram in the spin- 
sector is indicated in Fig. 2. 
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der phase transition. In the BCS-like regime, however, 
the nature of the transitions or sequence of transitions 
between these states is less clear. 

Rather than attempting to fully characterize this evo- 
lutionary process, we will content ourselves instead with 
determining the spin Lagrangian of the d-wave states 
themselves. Following the logic of the previous subsec- 
tion, consider the simplest modification of Eq. ^9| appro- 
priate for square anisotropy, 



+d 2 (\D 2 $ a \ 2 + |£>|$ CT | 2 ) 
+C cs {a a )-C CS {A°). 



di|D 2 $ 



(41) 



Here d 2 is a measure of the anisotropy, and low-energy 
excitations are pushed to finite momentum if c becomes 
negative. For c < and d 2 < 0, the lowest energy modes 
occur at four points k = (±K, 0), (0, ±K), precisely as 
expected for the low-energy excitations in a d xy super- 
conductor (the momenta appropriate for d x 2_ y 2 are ob- 
tained for d2 > 0)1 We fix c < and d 2 < and consider 
decreasing r a . 

Once again, it is appropriate to focus on the four 
slowly- varying fields $ JS (j = 1, 2, s = ±), defined by 



SK; 



(42) 



w here Ki = ( K,0), K 2 = (0,K), and K = 
\J\c\/ (2(di + d 2 )). As for the p-wave case, an effective 
theory can be developed for the $j S fields. It takes the 
form 

£e// = £ lA^-sl 2 + v 2 F \D^ JS \ 2 + v 2 A \e or D 3 ,^ os \ 2 



Fig. 2: The phase diagram as a function of increas- 
ing lattice potential with uniaxial anisotropy, Vi a ttico , and 
short-range attractive pairing force, ^attractive- 



B. square symmetry and d- wave 

We have seen how a BCS-like p x state with gapless 
modes at non-zero momentum can emerge from the spin 
boson Lagrangian in Eq. |2^ in the presence of uniaxial 
anisotropy. If the underlying crystal lattice has square 
symmetry, however, the p x or p y states are much less 
likely, and we expect instead that d-wave pairing is fa- 
vored. Indeed, two particles in a fourfold rotationally 
invariant lattice potential interacting with a hard core 
will generally have a d-wave ground state in the limit of 
a strong potential. For a weak anisotropy, however, the 
p x ± ip y states are likely lower in energy. One thereby 
expects a transition upon increasing lattice coupling be- 
tween p-wave and d-wave pairing. In the molecular limit 
this presumably occurs as a level crossing, i.e. a first or- 



(43) 



Neither the precise values of vf, «a nor the form of W 
is critical to this discussion. What is significant, how- 
ever, is the fact that Eq. ^ takes the form of four in- 
teracting relativistic complex bosons coupled to a sin- 
gle U(l) Chcrn-Simons gauge field. As before, the ef- 
fect of this gauge coupling is to attach identical flux 
to all spin quanta, transmuting these bosons into four 
species of Dirac fermions. Furthermore, the spin Hall 
conductivity a xy = 1 in the massive phase (i.e. for 
r a — c 2 /(4(g?i + d 2 j) > in mean- field theory), as de- 
termined by the last term in Eq. |^|. This is in agree- 
ment with the value in the p x ± ip y phase. Here, this 
value gives some indication of the structure of the signs 
of the four Dirac masses in the refermionized version of 



Eq. 43 , Each massive Dirac equation gives a contribution 



of ±1/2 to <J xyl depending upon the sign of its mass term, 
so apparently there must be three positive and one neg- 
ative (or vice versa) masses in the Dirac theory. As the 
quadratic term in W is tuned to zero, all four bose (and 
hence Dirac) fields are expected to go critical. This is 
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the natural candidate for a time-reversal invariant point, 
and we speculate that it describes the nodal quasiparti- 
cles (nodons) in a d~y. superconductor and relatives such 
as the nodal liquid.tll The appropriate fermionic repre- 
sentation for the nodons in these states was derived in- 
dependently in Ref. [H]. For d xy symmetry, it can be 
obtained as before from Eq. ^6| by defining 

c i( x ) = E d M elsKj " x ' ( 44 ) 

The (massless) nodon Lagrangian is then 

Cnodon = E + sv F r z idi + svAT x id 2 ] d ls 

+d\ s [d T + sv F T z id 2 + sv&T x di] d 2s , (46) 

where as in earlier equations the f matrices act in the 
a =T,| subspace. For an alternate, explicitly SU(2) 
invariant formulation, see Ref. |ll|. Eq. ^ has the de- 
sired form of four Dirac equations (for s = ±, j = 1,2). 
The identification of Eq. is supported by the excellent 
correspondence between the quantum numbers and mo- 
menta of the gapless modes of Eq. ^ and those of the 
nodons. In either case, there are four conserved U(l) cur- 
rents, the charges (time components) of which are chiral 
spin densities, i.e. spin densities for particles with mo- 
menta along ±Ki, 1LK2. Mass terms taking Eq. [l6] away 
from criticality can also be added, and take the form of 
d\ s T v d is operators. The explicitly break time-reversal in- 
variance, so that as for the p x phase, we expect the d xy 
(and analogously d x 2_ y 2) state to be locally stable. 

Given the complexity of the arguments in this sec- 
tion, it seems appropriate to summarize what has been 
learned. We have studied how momentum space struc- 
ture emerges from the bosonic Ginzburg-Landau theory 
of the spin sector. In doing so, we have not assumed (as 
in previous work on the nodal liquid) local superconduc- 
tivity, but proceeded instead on very general grounds. 
Once the soft-modes of $ CT move to non-zero momenta, 
any incipient critical points can invariably be expressed 
in terms of multiple Dirac fields. When the associated 
Dirac masses vanish, a time-reversal invariant lagrangian 
is possible, and two such theories were identified with 
the p x and d xy nodal states. Furthermore, uniaxial and 
square lattice anisotropies were seen to favor appropriate 
critical states, even in very naive treatments of the Lan- 
dau theory. These arguments provide a partial deriva- 
tion of Eq. ^6] for a generic time-reversal invariant d-wave 
(superconducting or nodal liquid) state, contingent only 
on the original hard-core assumption used to allow the 
Chern-Simons flux attachment. 



VI. DISCUSSION 

In this paper, we have suggested that spin-charge sep- 
aration is a generic consequence of strong repulsion be- 
tween electrons in two-dimensions. We are driven in- 
exorably to this conclusion by the following logic. (1) 
We note that if there are repulsive interactions which are 
strong enough to prevent electron trajectories from inter- 
secting, then we may transmute the electrons into hard- 
core bosons interacting with a Chern-Simons gauge field 
which attaches flux to spin. (2) Since the up- and down- 
spin boson currents are separately conserved, they can be 
written as the curls of two auxiliary gauge fields. (3) Us- 
ing (1) and (2), we formulate an equivalent dual theory 
which is of Ginzburg-Landau form. The auxiliary gauge 
fields are minimally coupled to vortex fields. (4) The 
Ginzburg-Landau theory contains a Z 2 symmetry which, 
if broken by vortex condensation, leads to spin-charge 
confinement, translational symmetry-breaking, and 'con- 
ventional' ordered phases such as the AF and CDW. (5) 
In order to study phases with spin and charge physics 
at different scales and without translational symmetry- 
breaking, we contemplate ^-symmetric phases in which 
vortex pairs condense. A phase diagram of spin-charge 
separated states is the upshot of crossing the rubicon fed 
by these five tributaries. 

By eschewing a conventional momentum-space ap- 
proach which assumes a Fermi surface, we have con- 
structed an effective field theory which does not fall un- 
der the rubric of Fermi liquid theory. The basic ex- 
citations of oui'|-theory are topological solitons in vor- 
tex condensatescJ. They are in no sense adiabatically 
connected to the electron and hole excitations of a free 
Fermi gas. They are also rather different from the 'holon' 
and 'spinon' conceptSj-which are introduced to solve the 
Gutzwiller constraint^. These objects are strongly- 
coupled and do not appear to be soliton-like in charac- 
ter. The topological spinless charge e excitation in our 
Ginzburg-Landau theory (a "holon" ) has an antiparticle 
with opposite charge, and at low energies can decouple 
from the neutral spin 1/2 Fermionic "nodon" excitation. 
Moreover, the "holon" already exists as a finite energy 
excitation within the Mott insulator - doping is not re- 
quired. 

Contrasting the present work with our earlier construc- 
tion of the nodal liquid, we see that the "holon" field is 
identical to the charged soliton of the nodal liquid but 
the nodon - which descended from an assumed quasipar- 
ticle at the nodes of a d x 2_ y i superconductor - is now 
a concept which naturally arises from the spatially non- 
uniform softening of a vortex-anti-vortex pair field. Our 
nodon and holon fields are properly seen as analogous to 
the charge and spin solitons of the one-dimensional elec- 
tron gas or the fractionally-charge quasiparticles in the 
quantum Hall effect. 

Despite the surprising ease with which our dual 
Ginzburg-Landau formulation captures spin-charge sep- 
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aration and superconductivity, the Fermi liquid phase 
seems to be missing. Generally, a dual vortex descrip- 
tion of a Fermi liquid is possible, as illustrated nicely 
for the case of spinless electrons. After transforming to 
spinless bosons via Chern-Simons and implementing a 
duality transformation, one readily obtains a simple dual 
Ginzburg-Landau theory This theory closely resembles 
Eq. |^, but with only a single vortex field, <E>, which is 
minimally coupled to a single gauge field with Chern- 
Simons dynamics. But more importantly, the theory is 
non-relativistic (i.e. there is a non-zero chemical poten- 
tial) and the gauge field necessarily has a non-zero aver- 
age, as the dual flux equals the number of electrons. So 
a vortex "vacuum" phase - the p x + ip y superconductor 
for spinful electrons - is not accessible without spin. The 
Wigner crystal phase of spinless electrons corresponds 
simply to condensing the single vortex field, ($) 0. In 
the Fermi liquid phase this vortex field must remain un- 
condensed, but with the vortices in a fluid state. This 
fluid of vortices presumably coexists with the fluid of 
particles (the Chern-Simons bosons which are the dual 
flux tubes) - the particle motion acting to scramble the 
vortex phase and vice versa. 

By analogy, a description of a spinful Fermi liquid via 
the dual Ginzburg-Landau theory (Eq. ||) presumably re- 
quires an uncondensed but "critical" fluid of both up and 
down spin vortices. But this analogy is problematic for 
two reasons: Firstly, both vortex fields, <I> Q , are rela- 
tivist^ with an equal number of positive and negative 
circulation vortices. But more importantly, in Eq. |8| the 
Chern-Simons term couples to the electron spin, effec- 
tively mediating a long-range statistical interaction be- 
tween spin up and spin down electrons. In a spinful 
Fermi liquid, the up and down spin quasiparticle exci- 
tations should essentially propagate independently. It is 
for these reasons that we suspect the impossibility of de- 
scribing a Fermi liquid phase in the dual formulation. 

Why should the dual Ginzburg-Landau theory pre- 
clude a Fermi liquid, given that the only explicit assump- 
tion needed for its derivation was that opposite spin elec- 
tron coordinates never coincide? In principle, the com- 
posite boson theory in Eq. || is exactly equivalent to a 
hard-core interacting fermion model, and therefore has 
a Fermi liquid phase in the presence of hard-core inter- 
actions only. The latter conclusion can be understood 
from the two-body problem of a small hard-core scat- 
terer, from which it is seen that the s-wave scattering 
amplitude is small when the hard-core radius a <C hp. 
The subsequent manipulations leading to the dual for- 
mulation, however, do not treat this constraint carefully. 
Since the flux attachment procedure itself adds relative 
angular momentum to the wavefunction (c.f. Eq. |l[), to 
accomodate a significant s-wave component requires in- 
corporating p-wave effects in the composite bosons. Such 
non-local physics is easily missed, and presumably re- 
quires a much more subtle treatment of the duality trans- 
formation condition to capture it. When the repulsive 
interactions are sufficiently strong, however, we suspect 



the dual lagrangian to be adequate. 

As we have argued, the presence of electrons spinning 
around one another is tantamount to significant (finite 
angular momentum) pairing correlations. The kinetic en- 
ergy clearly favors lower angular momentum, which sug- 
gests the predominance of p x +ip y pairing, at least in the 
absence of significant ionic potentials. Naively, this rea- 
soning might suggest that quantum Hall systems with 
spinful electrons could exhibit high-temperature super- 
conductivity. However, the presence of the strong orbital 
magnetic field presumably precludes this possibility. But 
there are quantum Hall systems which apparently do ex- 
hibit high temperature "pseudo-spin" superfluiditjcS. In 
particular, double-layer quantum Hall systems with total 
filling v = 1 do exhibit evidence of a "transverse super- 
fluid" phase, with superfluid currents in the two layers 
flowing readily in opposite directions (effectively negat- 
ing the effects of the magnetic field). It has been sug- 
gested that this phase will disorder via a finite tempera- 
ture Kosterlitz-Thouless transition, with Tkt in the 1/2 
Kelvin range - a superfluid transition driven by Coulomb 
repulsion. 

The rather close analogy between this quantum Hall 
system and the 2d p x + ip y superconductor (under the 
exchange of spin and charge) suggests a promising route 
for studying energetics. In the quantum Hall effect, vari- 
ational wavefunctions have provided a powerful means of 
comparing various candidate ground states. Interchang- 
ing charge and spin, we therefore propose the following 
candidate wavefunction for the p x + ip y superconductor: 

*uT = n> - - «i)n^* -%) e " 7Ei " |a+iu " |a - 

i>3 i,3 

(47) 

Here, z% = x% + iyi denotes the complex coordinate of an 
up spin electron, and Wi are down spin electron coordi- 
nates. The overbar denotes complex conjugation. As for 
quantum Hall wavefunctions the electron density is set 
by the exponential terms (with 7 proportional to the 2d 
density). But this wave function involves both z and it's 
complex conjugate z. In this way, electrons moving with 
large angular momentum are avoided. In a slight abuse 
of quantum Hall notation, we refer to the wavefunction 
in Eq. |47] as a "(1,1,1)" state. Readers familiar with 
Chern-Simons theory will recognize that Eq. ^ encapsu- 
lates the universal content ( "K- matrix" ) of the effective 
field theory in Eq. [l3|, and is therefore a faithful repre- 
sentation of the p x + ip y state. To make the "pairing" 
more explicit, it is instructive to use the Cauchy identity 
to re-express this wavefunction as 

* UT = Dot f-J— ) HI* - ^| 2 e-5>l 2 +W 2 . 

(48) 

Eq. ^ displays as the product of a BCS wavefunc- 
tion (projected onto a state of definite electron number) 
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with pair wavefunction 1/z and a Jastrow factor. The 
Jastrow factor keeps up and down spin electrons apart 
appropriate for repulsively interacting 2d electrons, but 
the first term evidently encapsulates pairing correlations 
adequate for superconductivity. It will be interesting to 
explore the energetics of V I' 11 - in realistic models of re- 
pulsively interacting Fermions. 



T 




x 

Fig. 3: Illustrative and schematic phase diagram of 
cuprate superconductors considered within the present 
dual Ginzburg-Landau framework. The (unconventional) 
Antiferromagnet, Holon Lattice, and SuperConducting 
phases are indicated by AF*, HL, and SC, respectively. 

In the canonical approach to strongly interacting elec- 
tron systems near a Mott insulating phase, the very first 
step is to project onto a simpler tight binding model, of- 
tentimes with one orbital per unit cell. When modeling 
the undoped cuprate superconductors a further projec- 
tion to a reduced Hilbert space with one electron per site 
is usually adopted. The resulting spin Hamiltonian is 
much more tractable than the full system of interacting 
electrons, but we maintain that very important physics is 
irretrievably lost under these projections. For instance, 
the spin-charge separation that we access readily within 
our dual Ginzburg-Landau formulation of strongly in- 
teracting electrons in the 2d continuum is certainly not 
present in the Heisenberg spin Hamiltonian. Put simply, 
it is exceedingly difficult to ascertain if spin and charge 
separate after projecting away the charge. We can, how- 
ever, within our continuum approach, still describe Mott 
insulating physics by including a commensurate ionic po- 
tential which locks the charge order. In this way it is pos- 
sible to describe Mott insulators which have gapped exci- 
tations and exhibit spin-charge separation. Indeed much, 
if not all, of the interesting physics accessible within the 



dual Ginzburg-Landau formulation is inaccessible, and 
probably not present, in the overworked t — J model. 

An apparent limitation of our approach is that we do 
not have a specific microscopic model which is described 
by our effective field theory and we do not know how 
strong the interactions must be to invalidate Fermi liq- 
uid theory. Nevertheless, we expect that, since many of 
the phases have excitation gaps, they should be stable 
to small perturbations and there should be a universality 
class of models which exhibit the same "universal" and 
robust properties (such as spin-charge separation). In- 
deed, as emphasized previously, the only formal require- 
ment for obtaining the dual formulation is the "hard- 
core" part of the repulsion. 

The physics of the Ginzburg-Landau formulation 
shares some tantalizing similarities with the cuprate 
high-T c materials. As discussed above, the crucial as- 
sumption of unbroken Z2 symmetry appears most natural 
for systems with a large native disparity between charge 
and spin ordering scales. This is indeed the case for 
the cuprates near their half-filled Mott-insulating states. 
Precisely at half-filling, charge fluctuations begin to be- 
come quantized (acquire a gap) at very high tempera- 
tures of order eV (the Hubbard U), while for spins the 
characteristic ordering energy scale is significantly lower 
of order J (with local moment formatio n s ec uring some- 
what higher). As in our previous work,! 1 HI 12 ! we hypoth- 
esize that this Mott phase comprises a spin-charge sep- 
arated insulator described by Eq. [l7]. As the electron 
density at half-filling is commensurate with the underly- 
ing CuC>2 plane periodicity, the charge sector (as well as 
the spin sector) is in effectively zero dual magnetic field. 

The Ginzburg-Landau theories for charge and spin 
suggest a behavior in the temperature-doping plane il- 
lustrated in Fig. 3. At half-filling, both vortex fields 
have zero external flux and make transitions from their 
"normal" states at high temperatures to their "Meiss- 
ner" states at low temperature. The associated mean- 
field transition temperatures, roughly T? MF ~ U and 
T° MF > J, are shown in Fig. 3. Below T° MF , the spin 
boson <& a begins to develop amplitude fluctuations, rep- 
resenting local moment formation. At somewhat lower 
temperatures this amplitude softens particularly near 
(±7r/2, ±7r/2), and the refermionized Dirac fields subse- 
quently order into an AF* phase (see Ref. |l2]for details). 
Doping x introduces a dual external flux V x a p oc fox 
into the charge sector only - the spin boson $ CT is largely 
unaffected and in particular T° MF presumably decreases 
only weakly. The dual flux in the Ginzburg-Landau the- 
ory for $ p introduces a dual mean-field U H C 2 V line, or 
rapidly decreasing Tf! MF {x). Within a mean-field treat- 
ment a holon lattice phase would be expected below this 
line, in direct analogy with the Abrikosov flux lattice. 
But with fluctuations the holon lattice phase should be 
separated from T£ (x) by a crossover regime analo- 
gous to the strongly-fluctuating "vortex liquid" state in 
type II superconductors. In the cuprate context, this is a 
regime of strong dynamical charge fluctuations and can 
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be thought of as a "holon liquid" , comprised of charge 
e bosons. The T[! (x) line then represents a crossover 
from a metallic phase above (with "unquantized" charge) 
to the holon liquid which manifests (dynamical) charge 
"quantization" in units of e (cf. to dynamical "flux 
quantization" in the vortex liquid). Remarkably, our 
Ginzburg-Landau formulation suggests that spin-charge 
separation, effectively present below T ( f MF (x), can occur 
on very high interaction energy scales (eg. of order U). 
Upon further cooling the holon liquid one expects these 
bosons to condense, provided their density is sufficiently 
incommensurate with the underlying crystal potential to 
avoid charge ordering (into, for example, a holon lattice 
phase). This is the superconducting state, expected to 
be d-wave with a strong four-fold ionic potential. The 
predominant effect upon cooling through T° MF above 
the superconducting phase, should be a reduction of low 
energy spin fluctuations and nodal formation in the elec- 
tron spectral function, with a lesser effect in the charge 
sector due to weak spin-charge couplings. Although we 
emphasize that this is very much a preliminary applica- 
tion of the ideas of this paper, the picture in Fig. 3 is 
suggestive. 

Our primary conclusion concerning the ubiquity of 
spin-charge separation and superconductivity driven by 
very strong repulsion has potential implications for a 
much broader class of other strongly interacting systems. 
Besides the cuprates, other systems include: the heavy 
fermion superconductorsu, quasi-one dimensional organic 
superconductors, low carrier 2DEG's with very laptte.r s in 
semiconductor MOSFET's and heterostructuresQca, su- 
perconductivity, in Sr2RuOi with possible p x + ip v pair- 
ing symmctryE3, the normal and superfluid phases of_3- 
He (the A phase with a p x + ip y pairing symmetryo) 
and perhaes-,most intriguingly the magnetic states of 
solid 3-HeO'E2l. In many of these systems one is also 
very much interested in the full three dimensional limit, 
particularly for 3-Hc. Unfortunately, the Chern-Simons 
approach transforming fermions into bosons by flux at- 
tachment is restricted to strictly two-dimensional sys- 
tems. But it is possible to transform between fermions 
and bosons in three dimensions by binding "statistical" 
magnetic monopoles to the particlesEJ. Unfortunately, 
this introduces an unphysical internal statistical mag- 
netic field (in contrast to the pure gauge coupling within 
Chern-Simons theory). But by attaching monopoles to 
spin, the monopole fields from the up spin electrons and 
the antimonopole fields from the down spin electrons 
will largely cancel (exactly on the average). Moreover, 
particle-vortex duality transformations are also possible 
in three-dimensions (ie. "electric- magnetic" duality), so 
it should be possible to obtain an entirely bosonic (but 
approximate) dual description of 3d electrons with a 
'non-crossing' constraint. Perhaps this approach might 
be useful in modeling some 3d strongly correlated sys- 
tems. 

If, as we have suggested, strongly-interacting spin-1/2 
fermions do not form a Fermi liquid, then our effec- 



tive field theory represents a new paradigm for corre- 
lated electron behavior. If, as we have further hypoth- 
esized, superconductivity is a prevalent attribute of the 
phases ensconced in our theory, then it is a paradigm 
which includes a new route to superconductivity. For 
these reasons, we submit that that our scenario could 
have far-reaching implications for the cuprate supercon- 
ductors and other strongly- interacting electron systems. 
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